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Abstract 

Inspired by the possibility of high-scale super symmetry breaking in the string landscape where 
the cosmological constant problem and the gauge hierarchy problem can be solved while the strong 
CP problem is still a challenge for naturalness, we propose a super symmetric KSVZ axion model 
with an approximate universal intermediate-scale (~ 10 11 GeV) supersymmetry breaking. To pro- 
tect the global Peccei-Quinn (PQ) symmetry against quantum gravitational violation, we consider 
the gauged discrete Zn PQ symmetry. In our model the axion can be a cold dark matter candidate, 
and the intermediate supersymmetry breaking scale is directly related to the PQ symmetry break- 
ing scale. Gauge coupling unification can be achieved at about 2.7 x 10 16 GeV. The Higes mass 
range is 130 GeV to 160 GeV. We briefly discuss other axion models with high-scale supersymmetry 
breaking where the stabilization of the axion solution is similar. 
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I. INTRODUCTION 



There exists a great variety of structures for the superstring/M theory vacua, the ensemble 
of which is called the stringy "landscape" With the "weak anthropic principle" ^, 
this proposal may not only give the first concrete explanation of the very tiny value of the 
cosmological constant which can only take discrete values, but also solve the gauge hierarchy 
problem. In particular, the supersymmetry breaking scale can be high if there exist many 
supersymmetry breaking parameters or many hidden sectors ^, ^ | . Although there is no 
definite conclusion that the string landscape predicts high-scale or TeV-scale supersymmetry 
breaking e 3], it is interesting to study models with high-scale supersymmetry breaking 

If supersymmetry is indeed broken at a high scale, the breaking scale can range from 1 
TeV to the string scale. For simplicity, we consider three representative cases: (1) string- 
scale supersymmetry breaking, (2) intermediate-scale supersymmetry breaking, and (3) TeV- 
scale supersymmetry breaking. The TeV-scale supersymmetry has been studied extensively 
during the last two decades, so we do not discuss it here. 

Although the original motivation for the New Minimal Standard Model (NMSM) is 
different from the above string landscape argument, the NMSM provides a concrete model 
with string-scale supersymmetry breaking. With only six extra degrees of freedom beyond 
the minimal Standard Model (SM), the NMSM incorporates the new discoveries of physics 
beyond the minimal Standard Model: dark energy, non-baryonic dark matter, neutrino 
masses, as well as baryon asymmetry and cosmic inflation, by adopting the principle of 
minimal particle content and the most general renormalizable Lagrangian. The NMSM is 
free from excessive flavor-changing effects, CP violation, too-rapidproton decay, problems 
with electroweak precision data, and unwanted cosmological relics . Moreover, gauge cou- 
pling unification can be achieved by considering threshold corrections |6[ , or by introducing 
suitable additional particles, similar to the axion models with string-scale supersymmetry 
breaking that will be discussed in Section V. The unitarity constraints on couplings are 
weaker than the stability and triviality constraints [2J. 

For high-scale supersymmetry breaking, Arkani-Hamed and Dimopoulos proposed the 
split supersymmetry scenario where the scalars (squarks, sleptons and one combination of 
the scalar Higgs doublets) have masses at an intermediate scale, while the fermions (gauginos 
and Higgsinos) and the other combination of the scalar Higgs doublets are still at the TeV 
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scale |4]. Gauge coupling unification is preserved and the lightest neutralino can still be 
a dark matter candidate. In addition, most of the problems with supersymmetric models, 
for example, flavor and CP violations, dimension-5 fast proton decay and the stringent 
constraints on the lightest Higgs mass, are solved. The consequences of split supersymmetry 
have been studied in Refs. a, |9(. 

On the other hand, unlike the cosmological constant problem and the gauge hierarchy 
problem, the strong CP problem is still a concern in the string landscape lQj . In the 
Standard Model, the 9 parameter is a dimensionless coupling constant which is infinitely 
renormalized by radiative corrections. There is no theoretical reason for 9 to be smaller 
than 10 -9 as required by the experimental bound on the electric dipole moment of the 



neutron 
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There is also no known anthropic constraint on the value of 9, i.e., 9 may 
be a random variable with a roughly uniform distribution in the string landscape [l0|. An 
elegant and popular solution to the strong CP problem is provided by the Peccei-Quinn 
(PQ) mechanism Q], in which a global axial symmetry U(1)pq is introduced and broken 
spontaneously at some high energy scale. Then 9 is promoted to a dynamical field, 9 = a/f a , 
with an effective potential for this field induced by non-perturbative QCD effects. The axion 
a is a pseudo-Goldstone boson from the spontaneous U(1)pq symmetry breaking, with a 
decay constant f a (of similar scale as the U(1)pq symmetry breaking scale). Minimization 
of the axion potential fixes the vacuum expectation value (VEV) of a, equivalently forces 9 
to be zero, thus naturally solving the strong CP problem. 

The original Weinberg- Wilczek axion |l4j is excluded by experiment, in particular by the 
non-observation of the rare decay K — > ti + a ^j] . There are two viable "invisible" axion 
models in which the experimental bounds can be evaded: (1) the Kim-Shifman-Vainshtein- 
Zakharov (KSVZ) axion model, which introduces a SM singlet and a pair of extra vector-like 
quarks that carry £ 7(1 ) pq charges while the SM fermions and Higgs fields are neutral under 
U(1)pq symmetry (2) the Dine-Fischler-Srednicki-Zhitnitskii (DFSZ) axion model, in 

which a SM singlet and one pair of Higgs doublets are introduced, and the SM fermions and 

I I 

Higgs fields are charged under U(1)pq symmetry [lfil ]. 

From laboratory, astrophysics, and cosmology constraints, the U(1)pq symmetry breaking 
scale f a is limited to the range 10 10 GeV < f a < 10 12 GeV [11]. Light axions can be 
produced in stars, and part of the star energy can be carried away by these axions. To avoid 
an unacceptable star energy loss, a lower bound is obtained on the axion decay constant, 
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fa > 10 10 GeV. During the course of cosmological evolution of the universe, the axions 
decouple early and begin to oscillate coherently. If f a is larger than about 10 12 GeV, at 
some point in the evolution the energy density in the coherent axion oscillations could 
exceed the critical energy density and over-close the universe. The invisible axion can be a 
good cold dark matter candidate 11]. 

In this paper we consider the axion models with high-scale supersymmetry breaking. Un- 
like split supersymmetry, we consider an approximately universal supersymmetry breaking, 
i.e., the mass differences among the squarks, sleptons, gauginos, Higgsinos and the A-terms 
are within one order of magnitude (within one-loop suppressions), which seems to be quite 



natural from flux- induced supersymmetry breaking [17 1. 

We construct a supersymmetric KSVZ axion model where we introduce a SM singlet S 
and two pairs of SM vector-like particles (Qx, Qx) an d (D x , D x )- The supersymmetry 
breaking scale and the Peccei-Quinn symmetry breaking scale are both around 10 11 GeV, 
and the axion can be a cold dark matter candidate. We give the Higgs potential and the SM 
Yukawa couplings below the intermediate scale, the superpotential above the intermediate 
scale, and the matching conditions for the couplings. In addition, we consider the discrete 
Zx Peccei-Quinn symmetry, which can be embedded into an anomalous U(1)a gauge sym- 
metry in string constructions where the anomalies can be cancelled by the Green-Schwarz 
mechanism 18]. This Zx discrete symmetry cannot be violated by the quantum gravita- 
tional interaction. In order that the contributions to the 9 term from the non-renormalizable 
operators can be less than 10~ 9 , we must require N > 10. We also show that a suitable 
quartic coupling for the singlet S can be obtained at one-loop due to its large Yukawa cou- 
plings to the extra vector-like particles. Thus, the intermediate supersymmetry breaking 
scale is directly connected to the Peccei-Quinn symmetry breaking scale. Moreover, us- 
ing two-loop renormalization group equation runnings for the gauge couplings and one-loop 
renormalization group equation runnings for the Yukawa couplings, we find that gauge cou- 
pling unification can be achieved at about 2.7 x 10 16 GeV and the Higgs mass range is from 
130 GeV to 160 GeV. Also, we calculate the Yukawa couplings for the third family of the 
SM fermions at the Grand Unified Theory (GUT) scale and comment on how to achieve 
Yukawa coupling unification. In our model, similar to the NMSM and split supersymmetry, 
most of the usual problems in the supersymmetric models are solved. If the i?-parity is an 
exact symmetry, the lightest supersymmetric particle (LSP) neutralinos can be produced 
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only non-thermally in a suitable amount, and their annihilations may have cosmo logical 
consequences. If i?-parity is broken, the tiny masses and bilarge mixings for active neutri- 
nos can be realized naturally at one-loop although we have to suppress the contributions to 
neutrino masses from the tree-level /^Ljif u terms in the superpotential. 

We also briefly discuss other axion models: the DFSZ axion model with intermediate-scale 
supersymmetry breaking where the stabilization of the axion solution and the Peccei-Quinn 
symmetry breaking are similar, the axion models with string-scale supersymmetry breaking, 
and the axion models with split supersymmetry. 

II. THE KSVZ AXION MODEL WITH INTERMEDIATE-SCALE 
SUPERSYMMETRY BREAKING 

We first specify our conventions. We denote the left-handed quark doublets, the right- 
handed up-type quarks, the right-handed down-type quarks, the left-handed lepton doublets, 
the right-handed leptons, and the Higgs doublet in the Standard Model as Ui, di, k, ti and 
H, respectively, where i is the family index. For the supersymmetric Standard Model, the 
SM fermions and Higgs fields are superfields belonging to chiral multiplets. We denote the 
left-handed quark doublets, the right-handed up-type quarks, the right-handed down-type 
quarks, the left-handed lepton doublets, the right-handed leptons, and one pair of Higgs 
doublets as Q iy U°, Df, L iy E%, H u and H d , respectively. The supersymmetry breaking scale 
is assumed to be at an intermediate scale in the range 10 10 GeV to 10 12 GeV. For simplicity, 
we assume that the gauginos, squarks, sleptons, Higgsinos, and one combination of the scalar 
Higgs doublets have a universal supersymmetry breaking soft mass M. 

We consider the KSVZ axion model jl^ with a discrete Peccei-Quinn symmetry Zx- 
This Zx symmetry can be realized after the anomalous U(1)a gauge symmetry is broken 
close to the string scale, where the U(1)a anomalies are cancelled by the Green-Schwarz 
mechanism j[18], as will be discussed in the next Section. In the KSVZ axion model, the 
SM fermions and Higgs fields are assumed to be neutral under the Zx PQ symmetry. We 
introduce one SM singlet S and two pairs of SM vector-like particles (Qx, Qx) an d {Dx, 
Dx)- Their quantum numbers under the SM SU(3)c X SU(2)l x U{1)y gauge symmetry 
and the Zx symmetry are specified in Table HI The axion potential from the QCD anomaly 
is generated when we couple the singlet S to these extra vector-like particles. In addition, 
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TABLE I: The quantum numbers for the superfields Qx, Qxi Dx, Dx, and S under the SU (3)c X 
SU(2)l x U(1)y x Zn symmetry. 



Superfields 


Quantum Numbers 


Superfields 


Quantum Numbers 


Qx 


(3, 2, 1/6, n Qx ) 


Qx 


(3, 2, -1/6, ngj 


D x 


(3, 1, -1/3, n Dx ) 




(3, 1, 1/3, n Sx ) 


S 


(1, 1, 0, n s ) 







the charges of Qx, Qx, Dx, D x , and S under the Z N symmetry satisfy the relations 

n Qx + % y +n s = mod N , 



n Dx + n-p + n s 



mod N . 



(1) 
(2) 



In our model the two pairs of the SM vector-like particles (Qx, Qx) an d (Dx, Dx) 
have masses comparable to the universal supersymmetry breaking soft mass M after the 
supersymmetry and Peccei-Quinn symmetry breakings. Below the scale M we assume that 
only one scalar Higgs doublet H is light as arranged by suitable fine-tuning. Thus, besides 
the gauge kinetic terms for all the fields, the most general renormalizable Lagrangian at 
tree-level is 

A 



C 



m 2 H H^H 



2! 



(iJtff) - [h^q jUl eH* + h'jyijdill + //;//,// + H.C. 



(3) 



where m 2 H is the squared Higgs mass, A is the Higgs quartic coupling, e = ia 2 where a 2 is 
the second Pauli matrix, and h^, hf-, and h\- are the Yukawa couplings. 

The Lagrangian in Eq. (jSJ) is obtained after the gauginos, squarks, sleptons, Higgsinos, and 
one combination of the scalar Higgs doublets from H u and in the supersymmetric model 
are decoupled at the scale M. Matching the Lagrangian in Eq. (jSJ) with the supersymmetric 
interaction terms for the Higgs doublets H u and Hd, we obtain the relevant Lagrangian in 
terms of the coupling constants in the SM at the scale M 



£susy = -f (Hla l H u + HyH d f - 9 f (hIH u - H\H C 



+ 



y%H u eUiqj - >/!jlI.i<'h'lj ~ y^Hdee^j + H.C 



(4) 



where a 1 are the Pauli matrices, #2 and gy are respectively the gauge couplings for the 
SU(2)l and U(l)y groups, and yfj, yfj, and y\j are the Yukawa couplings at and above scale 
M. Here we neglect the Yukawa couplings for the extra vector-like particles that will be 
given later in Eq. (J7|). 

From the scalar Higgs doublets in the chiral multiplets, H u and H^, we define the SM 
Higgs doublet combination as H = — cos f3tH^ + sin f3H u which is fine-tuned to have a 
small mass term, while the other combination (sinfleHj + cos (3H U ) has mass around M j^. 
Therefore, we obtain the coupling constants in Eq. (j3J) at the scale M from those in Eq. (0J) 
by the replacements H u — > sin (3H and Hd — > cos (3eH* 

[g 2 2 (M) + 3gj(M)/5\ 
A(M) = J= i cos 2 2/3, (5) 



hij(M) = (y%)*(M) sin/?, ^ e (M) = (^ e )*(M) cos/5 , (6) 

where we assume the SU(5) normalization for the gauge coupling of U(1)y, i-e., g\ = 
(5/3)^. 

Above the scale M the theory is supersymmetric. The superpotential in our model is 
W = vlH^UfQi - y%E d eD\Q 3 - y^eE^ + y Qx SQ x Q x 

+y Dx SD x D x + ysjjM-3 ' ( 7 ) 
m pi 

where yQ x , yo x an d ys are the Yukawa couplings for the extra particles, and Mpi is the 
Planck scale. The last term, ysS N j 'Mp[~ 3 ', is the lowest possible higher dimensional operator 
for the singlet S; this operator is related to the destabilization of the axion solution and 
suppressed by the Planck scale. Consequently, in the following discussions we do not consider 
the renormalization group evolution for the coupling ys and neglect it in our renormalization 
group equations. 



III. STABILIZATION OF THE AXION SOLUTION AND PECCEI QUINN 
SYMMETRY BREAKING 

Quantum gravitational effects, associated with black holes, worm holes, etc., are believed 
to violate all the global symmetries, while they respect all the gauge symmetries . These 
effects may destabilize the axion solution to the strong CP problem due to the violation 
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of the global Peccei-Quinn symmetry. However, after a gauge symmetry is spontaneously 
broken, there may exist a remnant discrete gauge symmetry which will not be violated by 
quantum gravity (^(J. Thus, a possible way to avoid the destabilization problem associated 
with quantum gravity is to identify the Peccei-Quinn symmetry as an approximate global 
symmetry arising from the broken gauge symmetry. 

In weakly coupled heterotic string model building, there generically exists an anomalous 
U(1)a gauge symmetry with its anomalies cancelled by the Green-Schwarz mechanism 18]. 
For Type II orientifold string model building, there may exist more than one anomalous 
U(1)a gauge symmetry whose anomalies can be cancelled by the generalized Green-Schwarz 
mechanism [21]. The anomalous U(1)a gauge symmetry is broken near the string scale 
when some scalar fields, which are charged under U(1)a, obtain VEVs and cancel the Fayet- 
Iliopoulos term of U(1)a- Then the D-flatness for U(1)a is preserved and the supersymmetry 
is unbroken [22j. Usually, there is an unbroken discrete Z N subgroup of the £/(1)a gauge 
symmetry, which is protected against quantum gravitational violation. We shall consider 
this Zn discrete symmetry as an approximate global U(1)pq symmetry [23]. 

For the gauge symmetry YliGi x U(1)a, the Green-Schwarz anomaly cancellation condi- 
tions from an effective theory point of view are 0, 13 

Aj A r , 



^-gravity 



Sgs , (8) 



ki 12 

where the Ai are anomaly coefficients associated with G\ x U(1)a, ki is the level of the 
corresponding Kac-Moody algebra, and Sqs is a constant which is not specified by low- 
energy theory alone. For a non-Abelian group, ki is a positive integer, while for the U(l) 
gauge symmetry, ki need not be an integer. All the other anomaly coefficients such as 
GiGjG k and [[/(l)^] 2 x d should vanish. 

In our model the gauge symmetry in which we are interested is SU(3)c X SU{2) L x 
U(l)y x U(1)a- Because the SM fermions and Higgs fields are not charged under U(1)a, 
there are no anomalies from them. Thus, we only need to consider the anomalies from the 
two pairs of SM vector-like fields (Qx, Qx) and (Dx, Dx) that involve at least one U(1)a- 

For simplicity, we assume that hq x = no x and n-q = n-^ x . Using the Georgi-Glashow 
SU(5) normalization j^, we find that the U(1)y[U(1)a] 2 anomaly vanishes and that the 
[U(1) Y ] 2 x U(1) A , [SU(2) l ] 2 x U(1) A and [SU(3) C ] 2 x U(1) A anomaly coefficients A u A 2 
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and As are 



5Al = A 2 = A 3 = -[n Qx +nQ x ) . (9) 

Thus, the anomalies can be cancelled by choosing the hi as follows 

h = i , A; 2 = fc 3 = 1 - (10) 

Usually, one only considers non-Abelian anomaly cancellations because the associated Kac- 
Moody level for £7(1) is not an integer in general, hence this condition is not very useful 
from a low-energy effective theory point of view |23L |2jJ . 

In addition, we can have models with k\ = k 2 = fc 3 = 1, although we have to introduce 
more particles. Let us give an explicit example. To ensure gauge coupling unification, we 
introduce two pairs of SM vector- like fields with quantum numbers ((3, 2, 1/6), (3, 2, —1/6)) 
and ((3, 1, —1/3), (3, 1, 1/3)). An alternative way to achieve gauge coupling unification 
is to introduce two chiral multiplets with quantum numbers (8, 1, 0) and (1, 3, 0), and one 
pair of SM vector-like particles with quantum numbers ((1, 2, —1/2), (1, 2, 1/2)). In 
these two cases, all the extra particles are neutral under U(1)a and have masses around 
10 11 GeV, and the gauge coupling unification scale is about 2 x 10 16 GeV. To realize the 
KSVZ axion and preserve the gauge coupling unification, we introduce at least one pair of 
SM vector-like fields with SU(5) quantum numbers, for instance, (5, 5) or (10, 10), that 
are charged under U(1)a and can couple to the singlet S via superpotential S55 or S'lOlO 
terms. Because these SM vector-like fields, which couple to S and give the axion potential 
from QCD anomaly, form a complete SU(5) representation, we have k\ = k 2 = k^ = 1. 
However, this model is not economical compared to above model because we need to add 
more particles. 

To ensure that the contributions to the 9 term from the non-renormalizable operators are 
less than 10~ 9 , we constrain the order N of the Zn symmetry. The most dangerous term is 
the following hard supersymmetry breaking term in the potential: 

S N 

V SB = Ay^s—j^ + E.C. , (11) 

M pl 

where A ys can be complex. This term would induce a non-zero contribution to 9 {27) of 

_ f N 
^ * AysVS - N \ A ■ (12) 
Mpi J\qcd 
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Here we assumed that CP violation is of order 1 so that there is no extra fine-tuning in our 
discussion. 

For numerical calculations, we first evaluate /„. Assuming that the ratio of the axion 
number density to the entro py d ensity has been constant since the axion acquired mass and 
started to oscillate, we have 12a 



l^rj • (13) 



f a y /6 ^200 MeV\ 3/4 

Using the dark matter relic density of Vt a h 2 ~ 0.1126 from the WMAP analysis 29] and 
taking A QCD = 200 MeV, we obtain f a ~ 1.54 x 10 11 GeV. 

Assuming that y$ = 1 and requiring that A6 < 10 -9 , we obtain 

N > 11, for A v , ~ 10 11 GeV 

VS , (14) 

N > 10, for A ys ~ 10 4 GeV 

* iriAmf ir\n r>h nco t~\v -f- Y\ o Tnorrnif nrlo A 

nsi 



where we take Mpi = 2.4 x 10 18 GeV. If the CP violation phase or the magnitude of A 



or the coupling ys are very small, the values A" can be smaller than these lower bounds. 

The above stabilization discussion is generic for supersymmetric axion models because the 
hard supersymmetry breaking terms can be induced from the non-renormalizable operators 
in the Kahler potential after supersymmetry breaking. For example, in our model we can 
have the following non-renormalizable operator in the Kahler potential 

8 = J d 4 xd*6 {^ ZSZ M Z N SN + H.C.) , (15) 

where yzs is a coupling constant, and Z is a spurion superfield which parametrizes the 
supersymmetry breaking via (Fz) ^ 0. For simplicity, we assume that the mass scale M* is 
equal to the Planck scale Mpi and (Fz)/Mpi is about 10 11 GeV. Therefore, we obtain 

A vs = V ~f - V ~f >< 10 4 GeV . (16) 

ys m pi ys 

Using y Z s ~ 1 and y s ~ 1, we obtain A^ > 10. 

In general, suppose we have several singlets Si that are charged under the Z^ symmetry 
whose VEVs are comparable to that of S, and suppose that the lowest dimensional operator 
allowed by the gauge symmetry and Zn symmetry is S no Y\ i S] l % where no is a positive 
integer and rij are non-negative integers. To ensure that the contribution to the 9 term from 
this operator is naturally less than 10 -9 , we require 

n +J2 n i ^ 11 or 10 • ( 17 ) 
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We next discuss Peccei-Quinn symmetry breaking. Because the soft mass and the VEV 
of the singlet S are about 10 11 GeV, we have to generate a suitable quartic coupling for S 
at one-loop from the two pairs of SM vector-like fields. This can be realized as long as dq x 
and ud x are of order 1. The one-loop effective potential for S is 

t r ~2| C |2 ^ 4 \q\±(\ ^Qxl^l 2 3\ 3 4 , c uf, Vd x \S\ 2 3\ 
V„, = -m s \S\ -^- 2 V Qx \S\ [l°t-Q^ - 2) - TW V °*W ( l0g -^ 2) 

^k^f^^^-l), , 18 ) 

where m| is the soft supersymmetry breaking mass for S, fnQ x is the soft mass for Q x and 
Q x , and rh 2 Dx is the soft mass for D x and D X ', Q is the renormalization scale, which is 
about f a here. Approximately, the quartic term for S is 

Veff D , (19) 

where 

For ?/q x = y Dx = 0.8 and m Qx = fn Dx = 5f a , we obtain X s = 0.17. 

We emphasize that for a supersymmetry breaking scale higher than 10 12 GeV, the VEV 
of S around f a ~ 1.54 x 10 11 GeV may not be natural. For a supersymmetry breaking 
scale lower than 10 10 GeV, one may not be able to generate the suitable quartic coupling 
for the singlet S naturally. Therefore, the Peccei-Quinn symmetry can be broken naturally 
if and only if the supersymmetry breaking scale is intermediate, i.e., from 10 10 GeV to 10 12 
GeV, and then the intermediate supersymmetry breaking scale is directly related to the 
Peccei-Quinn symmetry breaking scale. 

IV. GAUGE COUPLING UNIFICATION AND PHENOMENOLOGY 
A. Gauge Coupling Unification 



We next consider gauge coupling unification. The two-loop renormalization group equa- 
tions for the gauge couplings and the one-loop renormalization group equations for the 
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Yukawa couplings are given in Appendix A. The renormalization group equation runnings 
are logarithmic and we shall neglect the uncertainties from threshold corrections. For sim- 
plicity, we consider the universal supersymmetry breaking mass M for the superpartners 
of the SM particles, and we assume that the masses for the two pairs of SM vector-like 
superfields (Qx, Qx) an d (Dx, Dx) are the same and given by mx- 

In numerical calculations we use the following initial values at the Mz scale in Ref. j^jj]: 

a -\M z ) = 128.91 ±0.02, 

sin 2 9 W (M Z ) = 0.23120 ± 0.00015 , 

M z = 91.1876 ± 0.0021 GeV , 

v = 174.10 GeV, (21) 

along with a top quark pole mass m t = 178. ±4.3 GeV and a s (M z ) = 0.1182±0.0027 [32]. 
We choose m x = 1.0 x 10 11 GeV, and uq x = yr> x = 0.8 at the m x scale. The supersymmetry 
breaking scale M and the GUT scale M GUT are determined by requiring the gauge couplings 
to unify at M GUT . For the MS top quark Yukawa coupling, we use the one- loop corrected 
value 3J|, which is related to the top quark pole mass by 

m t = h t v(l + —JL - 2-^-) . (22) 

V 3 167T 2 167T 2 J K ' 

The unifications of the gauge couplings for tan f3 = 1.5 and tan /3 = 50 are shown in Fig.Q 
in which the supersymmetry breaking scale M is 1.4 x 10 10 GeV. The GUT scale Mqut is 
about 2.6 x 10 16 GeV and 2.7 x 10 16 GeV for tan/3 = 1.5 and tan/3 = 50, respectively. The 
supersymmetry breaking scale, the GUT scale, and the unified gauge coupling are almost 
independent of tan (5. 

We point out that if the gauginos and Higgsinos have masses of about M while the squarks 
and sleptons are about one order of magnitude heavier, the gauge coupling unification scale 
is the same as above because the squarks and sleptons form complete SU(5) representations. 

B. Higgs Mass 

To calculate the Higgs mass, we first evaluate the Higgs quartic coupling A at the scale 
M and then evolve it down to the Mz scale. The effective potential, with the top quark 
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FIG. 1: Gauge coupling unifications for tan/3 = 1.5 (left) and tan/3 = 50 (right). 



radiative correction included, is 



A 



eff 



m 



16vr 



2'"t 



log 



h 2 t (H^H) 3 
Q 2 2 



(23) 



We derive the Higgs mass by minimizing the effective potential. The scale Q is chosen to 
be at the Higgs mass. We show the predicted Higgs mass as a function of tan (3 in Fig. |21 
The upper, center and lower groups of curves correspond to the top quark mass being 
(178.0 + 4.3) GeV, 178.0 GeV and (178.0 - 4.3) GeV, respectively. Within each group, the 
upper and lower curves show the uncertainty due to the la uncertainty of a s . The lower 
curve is generated with a s = 0.1182 + 0.0027 and the upper with a s = 0.1182 — 0.0027. 

Below the scale M, the Higgs boson is just as in the SM. A SM Higgs boson in the mass 
range 130-160 GeV can be identified at the 3<r level of significance with 30 fb _1 integrated 



luminosity at the Tevatron 34J, 35]. At the LHC, the Higgs boson in this mass range can be 



discovered in the 77 and the r + r channels with 100 fb 1 of data 



C. Top, Bottom and Tau Yukawa Couplings at the GUT Scale 

Besides gauge coupling unification, it is also interesting to study whether there exists 
Yukawa coupling unification for the top (t) quark, the bottom (6) quark and the r lepton at 
the GUT scale. The ratios of the Yukawa couplings (yt/y T , Vt/Vb an d yb/yr) as a function 
of tan f3 are shown in Fig. El 
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FIG. 2: The Higgs mass versus tan /3. The group of solid curves are results for rat = 178.0 GeV, 
the dashed for rat = 178.0 + 4.3 GeV, and the dotted for rat = 178.0 — 4.3 GeV. Within each group, 
the center curve is generated with a s = 0.1182, the lower curve with a s = 0.1182 + 0.0027, and 
the upper with a s = 0.1182 - 0.0027. 



FIG. 3: The ratios of Yukawa couplings at the GUT scale versus tan/3. The solid curve is yt/yr, 
the dashed yt/yb, and the dotted yb/yr- 

From Fig. El we find that the t-T Yukawa couplings would be unified at tan/3 ~ 49, while 
the t—b or b-r Yukawa couplings cannot be unified. Thus, how to explain the t, b, and r 
Yukawa couplings at the GUT scale remains an interesting problem. This may be solved via 
flavor symmetry breaking or via particle mixings by introducing additional SM vector-like 
particles with similar SM quantum numbers. 
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D. Comments 



In our models, the gauge coupling unification scale is about 2.7x 10 16 GeV; therefore, there 
is no dimension-6 proton decay problem. Due to intermediate-scale supersymmetry breaking, 
we do not have the dimension-5 proton decay problem, the moduli problem, and the excessive 
supersymmetry flavor violation and CP violation, etc. In addition, if we want to explain 
neutrino masses via the see-saw mechanism 37 1 and baryon asymmetry via leptogenesis , 
we just need to introduce three right-handed neutrinos. Also, the SM fermion masses and 
mixings can be generated naturally via the Froggatt-Nielsen mechanism j^. 

If .R-parity is conserved in our model, the lightest supersymmetric particle, for example 
the neutralino, cannot be produced thermally, otherwise the universe will be over-closed. 
So we have to require that the reheating temperature be below 10 10 GeV. However, there 
is the interesting possibility that the LSP neutralino is produced non-thermally in suitable 
amount, and may still be a dark matter candidate. Then neutralino annihilation may have 
cosmological consequences. 

If i?-parity is violated, tiny neutrino masses can be generated naturally at one-loop. How- 
ever, for naturalness, we have to suppress the tree-level contributions to the neutrino masses 
from the fiiLiH u terms in the superpotential, where /ij must be smaller than 10 GeV j^j]. 
The neutrino masses and bilarge mixings deserve further detailed studies. 



V. THE OTHER AXION MODELS WITH HIGH-SCALE 
SUPERSYMMETRY BREAKING 

We now briefly discuss the other axion models with high-scale supersymmetry breaking, 
where the axion can be the dark matter candidate and the axion solution to the strong CP 
problem may be stabilized similarly to the discussions in Section UTTl However, the Peccei- 
Quinn symmetry can be broken naturally at about 10 11 GeV if and only if the supersymmetry 
breaking scale is intermediate. 



A. DFSZ Axion Model with Intermediate-Scale Supersymmetry Breaking 

Since we consider approximate universal supersymmetry breaking, i.e., the mass differ- 
ences among the squarks, sleptons, gauginos, Higgsinos and the A-terms are within one 
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order of magnitude (within one- loop suppressions), we must introduce additional particles 
to achieve gauge coupling unification. Similar to the discussions in Section III, we can in- 
troduce two pairs of SM vector-like particles with SM quantum numbers ((3, 2, 1/6), (3, 2, 
— 1/6)) and ((3, 1, —1/3), (3, 1, 1/3)). Alternatively, we can introduce two chiral mul- 
tiplets $8 an d $3 with quantum numbers (8,1,0) and (1,3,0), and one pair of the SM 
vector-like fields with quantum numbers ((1, 2, —1/2), (1, 2, 1/2)). The masses for these 
extra particles are around the intermediate scale 10 11 GeV, and gauge coupling unification 
can be achieved around 2.0 x 10 16 GeV. 

The stabilization of the axion solution is similar to that in Section III for the KSVZ 
axion model by introducing a discrete Peccei-Quinn symmetry and embedding it into 
the anomalous U(1)a gauge symmetry, whose anomalies are cancelled by the Green-Schwarz 
mechanism. As in Section III, in order to generate the quartic coupling for the singlet S, we 
need to couple S to the extra particles. For example, the superpotential for the additional 
particles in the second set of fields is 

W = y* 8 S® 8 <£> 8 + y$ 3 S®3®3 , (24) 

where the Yukawa couplings y$ s and y$ 3 are of order 1. Here $ 8 and $ 3 are charged under 
the Zn or U{1)a symmetry. 

We emphasize that in this model the intermediate supersymmetry breaking scale is still 
directly related to the Peccei-Quinn symmetry breaking scale. 



B. Axion Models with String-Scale Supersymmetry Breaking 



The axion models with string-scale supersymmetry breaking are the Standard Model plus 
the KSVZ or DFSZ axion. There are many ways to achieve gauge coupling unification by 
introducing various sets of particles with masses from the TeV scale to the GUT scale. For 
simplicity, we only consider the axion model with TeV-scale extra fermions whose masses 
can be protected by the chiral symmetry. We can introduce two pairs of SM vector-like 
fermions with quantum numbers (73, 2, 1/6), (3, 2, —1/6)) and ((3, 1, —1/3), (3, 1, 1/3)), 
similar to the extra particle content in Ref. 4l(. Gauge coupling unification can be achieved 



around 2 x 10 16 GeV 



4lj . For the KSVZ axion model with the second set of fields, we need 



to introduce SM vector-like particles which form complete SU(5) multiplets and couple to 
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the SM singlet S so that we can have the QCD anomalous Peccei-Quinn symmetry and 
preserve gauge coupling unification. 

The stabilization of the axion solution is similar to that in Section III. However, we need 
extra fine-tuning to keep the VEV of S around 10 11 GeV. 

C. Axion Models with Split Super symmetry 

For axion models with split supersymmetry, the stabilization of the axion solution and 
the Peccei-Quinn symmetry breaking are similar to those in Section III if the supersymmetry 
breaking scale is intermediate, i.e., from 10 10 GeV to 10 12 GeV. However, for axion models 
where the supersymmetry breaking scale is higher than 10 12 GeV or lower than 10 10 GeV, 
the stabilization of the axion solution is similar to that in Section III, while how to naturally 
break the Peccei-Quinn symmetry deserves further study. 

For the KSVZ axion model with intermediate-scale supersymmetry breaking, we can 
introduce one pair of vector-like fields with SU (5) quantum numbers 5 and 5, respectively, 
that can couple to the singlet S via the superpotential 555. As a result, gauge coupling 
unification is still preserved, and the stabilization of the axion solution and the Peccei-Quinn 
symmetry breaking are the same as those in Section III. For the DFSZ axion model with 
intermediate-scale supersymmetry breaking, we still need to introduce some adjoint or SM 
vector-like particles which form complete SU (5) multiplets that couple to the singlet S to 
generate its quartic coupling at one-loop. 

There are again two possibilities for the axion models with split supersymmetry: (1) im- 
parity is violated because the LSP neutralino need not be a dark matter candidate at all; 
(2) imparity is preserved, and both the LSP neutralino and the axion or dominantly one of 
them contributes to the dark matter density. 

VI. CONCLUSIONS 

There exists the possibility of high-scale supersymmetry breaking in string compactifi- 
cations where the cosmological constant problem and the gauge hierarchy problem can be 
solved. However, the strong CP problem may still be a serious complication for naturalness 
in the string landscape. Motivated by these considerations, we constructed a supersym- 
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metric KSVZ axion model where the supersymmetry breaking scale and the Peccei-Quinn 
symmetry breaking scale were both around 10 11 GeV and the axion could be a viable cold 
dark matter candidate. We considered the discrete Z N Peccei-Quinn symmetry that could 
be embedded into an anomalous U(1)a gauge symmetry in string constructions where the 
anomalies were cancelled by the Green-Schwarz mechanism. This Zj^ discrete symmetry 
could not be violated by quantum gravitational corrections. We found N > 10 is necessary 
to ensure that the contributions to the 9 term from the non-renormalizable operators are 
under control. We also showed that a reasonable quartic coupling for the singlet S could 
be generated from the one-loop corrections of extra vector-like particles. Thus, the interme- 
diate supersymmetry breaking scale was directly connected to the Peccei-Quinn symmetry 
breaking scale. In addition, using two-loop renormalization group equation runnings for 
the gauge couplings and one-loop renormalization group equation runnings for the Yukawa 
couplings, we showed that gauge coupling unification was achieved at about 2.7 x 10 16 GeV 
and that the Higgs mass was in the range 130 GeV to 160 GeV. We also calculated the 
Yukawa couplings for the third family of the SM fermions at the GUT scale and commented 
on the possibility of Yukawa coupling unification. In our model, due to intermediate-scale 
supersymmetry breaking, we did not have the dimension-5 operator induced proton decay 
problem, the moduli problem, and the problematic supersymmetry flavor and CP violations, 
etc. We also pointed out that if the -R-parity was an exact symmetry, the LSP neutralino 
could be produced only non-thermally and might have cosmological consequences. If the 
i?-parity was broken, the tiny masses and bilarge mixings for the left-handed neutrinos could 
be realized naturally at one-loop, although we had to suppress the tree-level contributions 
to the neutrino masses from the /ijLjif u terms in the superpotential. 

Furthermore, we briefly discussed the other axion models: the DFSZ axion model with 
intermediate-scale supersymmetry breaking where the stabilization of the axion solution 
and the Peccei-Quinn symmetry breaking were similar, the axion models with string-scale 
supersymmetry breaking, and the axion models with split supersymmetry. 
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APPENDIX A: RENORMALIZATION GROUP EQUATIONS 



In this Appendix, we give the renormalization group equations in the SM, and in our 
super symmetric KSVZ axion model. The gene ral formulae for the renormalization group 
equations in the SM are give n in Refs. 421 J^. and these for the supersymmetric models 



are given in Refs. 



44, 



45 



46|. 



First, we summarize the renormalization group equations in the SM. The two-loop renor- 
malization group equations for the gauge couplings are 

3 



T, B H9j- E d«Tr(h°th a 

3=1 



a=u,d,e 



(Al) 



where t = ln/i and \x is the renormalization scale. The g\, gi and g% are the gauge couplings 
for U(l)y, SU{2)l and SU(3)c, respectively, where we use the SU(5) normalization g\ = 
(5/3)g Y - The beta-function coefficients are 
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(A2) 



(A3) 



Since the contributions in Eq. ()A1|) from the Yukawa couplings arise from the two-loop 
diagrams, we only need Yukawa coupling evolution at the one-loop order. The one-loop 
renormalization group equations for Yukawa couplings are 



(47T) 



2 j f h u = h u dig] + - + A 2 ) 



(4tt) 2 ^ h d 
v ; dt 



(A4) 
(A5) 
(A6) 
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where 



I' I- 8 ! ' 



A 2 = Tr(3A ut A" + 3h d 'h d + h"h') . 

The one-loop renormalization group equation for the Higgs quartic coupling is 

d . _ _ - o /9 o „ o\ - 9/3 /i 2 



(A7) 
(A8) 



(47T) 2 ! A = 12A 2 - (|^ 2 + 9gl) A + \ (±g\ + \g\g\ + g() + 4A 2 A - 4A 4 , (A9) 



where 



Tr 



3{h u ^h u ) 2 + 3{h d ^h d ) 2 + {h e ^h e y 



(A10) 



Second, we summarize the renormalization group equations in our supersymmetric KSVZ 
axion model. The two-loop renormalization group equations for the gauge couplings are 

3 



(*0'1 ft = 9 3 ( a, + A6i) + _jL_ 



pB« + AB«) 9 J - £ d? Tr (j^V") 

^yL^Q* - d ? x ^*^J , (ah) 



where Abi and Af?^ are the contributions from the extra two pairs of the SM vector-like 
particles. The beta-function coefficients are 

'33 . \ . , /3 
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(A12) 

(A13) 

(A14) 
(A15) 



The one-loop renormalization group equations for Yukawa couplings are 



3»"V + !/V + 3Tr(uV) - £c? 9 . 

1=1 

y V + 3 y d V + Tr(3y d V + y^y e ) - £ c d g\ 



(A16) 
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where 
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